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CALIBRATION TECHNIQUES FOR IMAGING DEVICES 

FIELD 

The invention relates to color imaging and, more particularly, to calibration 
techniques for color imaging devices. 

BACKGROUND 

Calibration of an imaging device can significantly improve color accuracy of 
images rendered by the device. For example, imaging devices such as cathode ray tubes, 
liquid crystal displays, plasma displays, and various printing devices, are often calibrated 
to determine adjustments that can be applied to either color input or drive data applied to 
the device. In either case, the adjusted data can then be used to control the imaging 
device such that the ultimate rendition of the image has improved color accuracy. 
Calibration can be used to account for drift in the imaging device to improve color 
accuracy. 

For example, the calibration of a cathode ray tube (CRT) may involve attaching a 
measurement device, such as a colorimeter, to the display screen to measure color output 
of the CRT. The measured output can then be compared to analytical expected color 
values to determine the color errors. The determined errors can then be used to modify a 
lookup table (LUT) in a video card associated with a host computer so that input color 
data can be converted in a manner that adjusts for the determined errors. The 
effectiveness and accuracy of the calibration process can substantially impact color 
accuracy. 

Precise color accuracy is particularly important for color intensive applications 
such as soft proofing. Soft proofing refers to a proofing process that makes use of a 
display device rather than a printed hard copy. Traditionally, color proofing techniques 
have relied on hard copy proofing, where proofs are printed and inspected in order to 
ensure that the images and colors on the print media look visually correct. For instance, 
color characteristics can be adjusted and successive hard copy prints can be examined in 
a hard proofing process. After determining that a particular proof is acceptable, the color 
characteristics used to make the acceptable proof can be reused to mass-produce, e.g., on 
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1 5 calibrating an imaging device such as a display device, a printer or a scanner. The 
3} techniques may involve characterizing the imaging device with a device model, wherein 
yj an average error between an expected value of the device model and measured output of 
the image device is on the order of an expected error. The invention can achieve a 

balance between analytical behavior of the imaging device and measured output. In this 

UJ 

y, 20 manner, adjustments to image data may be more likely to improve color accuracy and 
less likely to overcompensate for errors that are expected. 

In various embodiments, the invention may be directed to methods of calibrating 
an imaging device. For example, a method may include characterizing the imaging 
device with a device model such that an average error between expected outputs 
25 determined from the device model and measured outputs of the imaging device is on the 
order of an expected error. The method may also include adjusting image rendering on 
the imaging device to achieve a target behavior. 

In another embodiment, the invention may be directed to a method that includes 
measuring outputs of the cathode ray tube for a subset of device values of the cathode ray 
30 tube, and choosing one or more parameter values of a device model, wherein the number 
of adjustable parameters is less than a number of measurements used to define the 
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a printing press, large quantities of print media that look visually equivalent to the 
acceptable proof. 

Soft proofing is desirable for many reasons. For instance, soft proofing can 
eliminate or reduce the need to print hard copies on media during the proofing process. 
5 Moreover, soft proofing may allow multiple proofing specialists to proof color images 
from remote locations simply by looking at display devices. With soft proofing, there is 
no need to print and deliver hard copy proofs to remote reviewers. Thus, soft proofing 
can be faster and more convenient than hard copy proofing. Moreover, soft proofing can 
reduce the cost of the proofing process. For these and other reasons, soft proofing is 
10 highly desirable. The ability to achieve precise calibration of soft proofing display 
devices is an important factor to achieving an effective soft proofing system. 

SUMMARY OF THE INVENTION 

In general, the invention is directed to various calibration techniques for 
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measured output of the cathode ray tube, and wherein an average error between expected 
outputs of the device model and the measured outputs is on the order of an expected 
error. The method may further include adjusting image data according to the device 
model to achieve a target behavior for the imaging device. 

In another embodiment the invention may be directed to a method that includes 
initializing a lookup table (LUT), adjusting settings of the cathode ray tube to 
substantially achieve a defined output, and measuring output for a number of color 
values. The method may also include choosing parameter values for a device model, 
wherein a number of adjustable parameters is less than a number of measured outputs, 
and generating entries for the LUT based on the device model. 

In another embodiment, the invention may implement a technique for biasing an 
output measurement by an amount sufficient to ensure that the output measurement is 
within a dynamic range of a measurement device. For example, a method may include 
measuring output of a display device, and displaying a substantially white trace during 
measurement to bias the output measurement. The trace may have a halo shape, or any 
other shape sufficient to properly bias the measurements. 

In other embodiments, the invention is directed to calibrated imaging devices or 
sets of calibrated imaging devices. For example, in accordance with the invention, a 
cathode ray tube, or a set of cathode ray tubes can be calibrated such that an average 
color error is approximately less than (0.75 delta e) from an analytical expected color 
output, and a maximum color error is approximately less than (1 .5 delta e) from the 
analytical expected color output. Furthermore, even more precise calibration, 
approaching a theoretical limits of analytic equations used to define device behavior can 
be achieved as described in greater detail below. 

Various aspects of the invention may be implemented in hardware, software, 
firmware, or any combination thereof. If implemented in software, the invention may be 
directed to a computer readable medium carrying program code, that when executed, 
performs one or more of the methods described herein. 

The invention is capable of providing a number of advantages. In particular, the 
invention can improve calibration of imaging devices. Moreover, improved calibration 
can facilitate the realization of color intensive applications such as soft proofing. In some 
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cases, the invention can be used to calibrate imaging devices such that measured errors of 
the imaging device are on the order of expected errors. For example, expected errors in 
the measurements may be caused by factors unrelated to the imaging device, such as 
errors introduced by the measuring device or the video card. The invention can achieve s 
5 balance between theory and measurement to ensure that adjustments to image data do not 
overcompensate for measured errors unrelated to the imaging device itself. In this 
manner, an imaging device can be calibrated such that an average color error is 
approximately less than (0.75 delta e) from an analytical expected color output. 
Additional details of these and other embodiments are set forth in the 
10 accompanying drawings and the description below. Other features, objects and 

advantages will become apparent from the description and drawings, and from the claims. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 is a front view of an exemplary imaging station including a cathode ray 

I s § 

15 tube calibrated according to the invention. 

BH FIG. 2 is a functional block diagram of an exemplary implementation of an 

SI . 

y3 imaging station according to the invention. 

j\ FIG. 3 is a flow diagram illustrating a calibration technique according to the 

iU invention. 

W 

|4 20 FIG. 4 illustrates one example of an imaging station displaying a white trace 

JEW. 

rf during calibration. 

FIGS. 5 and 6 are additional flow diagrams illustrating calibration techniques 
according to the invention. 

FIG. 7 is another functional block diagram of an exemplary implementation of an 
25 imaging station according to the invention. 

FIG. 8 is a block diagram of an exemplary soft proofing system including a 
number of imaging stations implementing the invention. 

DETAILED DESCRIPTION 

30 In the discussion that follows, many aspects of the invention are described with 

reference to the calibration of an imaging device in the form of a cathode ray tube (CRT). 
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However the invention is not necessarily limited in that respect. For example, techniques 
according to the principles of the invention may be readily applicable to other imaging 
devices, including other display devices such as liquid crystal displays, plasma displays, 
projection displays, and the like; printing devices such as printing presses, laser printers, 
ink-jet printers, dot-matrix printers, or any other printing device; and other imaging 
devices such as scanners. Accordingly, the detailed discussion is meant to be an 
exemplary description of one detailed embodiment in accordance with the invention. 

FIG. 1 is a front view of an exemplary imaging station 10. Imaging station 10 
includes an imaging device in the form of a cathode ray tube (CRT) 12. In addition, 
imaging station 10 includes a computer 14 that can receive image data and drive CRT 12 
according to the received image data. Imaging station 10 utilizes one or more of the 
calibration techniques described herein to improve rendering of color imagery on CRT 
12. In accordance with the invention, CRT 12 can be calibrated such that an average 
color error is approximately less than (0.75 delta e) from an analytical expected color 
output, and a maximum color error is approximately less than (1.5 delta e) from the 
analytical expected color output. More specifically, in some cases CRT 12 can be 
calibrated such that an average color error is approximately between (0.3 delta e) and 
(0.75 delta e) from the analytical expected color output, or even approximately between 
(0.3 delta e) and (0.4 delta e) from the analytical expected color output. In that case, the 
maximum color error may be approximately between (0.6 delta e) and (1.1 delta e) from 
the analytical expected color output, or even approximately between (0.6 delta e) and (0.8 
delta e) from the analytical expected color output. In other words, the invention can 
facilitate calibration accuracy of an imaging device that approaches a theoretical limit. 

Computer 14 may substantially conform to conventional computers used by 
graphic artists and other users in the creation of graphic imagery for electronic display or 
print production. For example, computer 14 may include a processor, memory, and an 
external storage device. A memory/bus controller and system bus typically couple the 
processor and memory, while one or more I/O controllers and an I/O bus couple the 
processor and memory to the storage device and CRT 12. Computer 14 may also include 
a user input device coupled to the processor an memory via an I/O bus. 
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The processor of computer 14 may take the form of a general purpose 
microprocessor and can be integrated with or form part of a PC, Macintosh, computer 
workstation, or the like. The user input device of computer 14 may include a 
conventional keyboard and pointing device such as a mouse, pen, or trackball, if desired. 
The memory of computer 14 may include random access memory (RAM) storing 
program code that is accessed and executed by the processor to carry out the calibration 
techniques described below. 

For example, program code implementing calibration techniques according to the 
invention can be loaded into the memory of computer 14 from the external storage device 
of computer 14, which may take the form of a fixed hard drive or removable media drive. 
The program code can be initially carried on computer-readable media such as magnetic, 
optical, magneto-optic or other disk or tape media. Alternatively, the program code may 
be loaded into memory from electronic computer-readable media such as electrically- 
erasable-programmable-read-only-memory (EEPROM), or downloaded over a network 
connection. If downloaded, the program code may be initially embedded in a carrier 
wave or otherwise transmitted on an electromagnetic signal. The program code may be 
embodied as a feature in an application program providing a wide range of imaging 
functionality. 

FIG. 2 is functional block diagram of an exemplary implementation of a imaging 
station 10 according to the invention. The various functional blocks may be implemented 
in hardware, or may be implemented in software which is executed in a processor within 
computer 14 as mentioned above. 

Imaging station 1 0 receives red-green-blue (RGB) image data as indicated at 
reference numeral 21 . Upon receiving RGB image data, the data can be manipulated by 
color matching module 22. In particular, color matching module may convert the RGB 
data using a color profile associated with the specific make and model of CRT 12. The 
converted data can then be sent through display driver 25 and video card 26 to ultimately 
drive the pixels of CRT 12 in a manner that yields a very accurate rendition of color 
images. 

To further improve color accuracy, imaging station 10 includes a calibration 
module 24 to calibrate CRT 12 in order to account for such things as drift in CRT 12. 
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Calibration module 24 can be invoked to measure output of CRT 12, and then, 
adjustments can be loaded into LUT 29 within video card 26. Video card 26 accesses 
LUT 29 to convert image values into drive values to be applied to CRT 12. These 
additional adjustments can be applied to the image data so that the ultimate rendition of 
CRT 12 has more accurate color. 

Calibration module 24 may invoke the calibration techniques described in greater 
detail below in order to achieve improved color accuracy. Typically, when a calibration 
technique is invoked, calibration module 24 prompts a user to attach a color measurement 
device to CRT 12. For example, the color measurement device is typically a high quality 
light detector that can be affixed to the display screen of CRT 12. Device number 
DTP92, commercially available from X-Rite Incorporated of Grand Rapids, Michigan is 
one suitable measurement device. After, the color measurement device has been affixed 
to the display screen of CRT 12, a user can initiate the calibration process. In that case, 
calibration module 24 performs a number of measurements on the output of CRT 12. 
After measuring the output of CRT 12, calibration module 24 can define a device model 
of CRT, and use the device model to generate adjusted values that can be loaded into 
LUT 29. In this manner, the calibration procedure can account for drift in CRT 12 and 
thereby improve color accuracy of images rendered on CRT 12. 

Calibration module 24 may begin with a presumption that CRT 12 is calibrated to 
a reference RGB color space having a perfect gamma curve behavior: 
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EQUATIONS 1, 2 A, 2B and 2C have normalized RGB values between 0.0 and 1.0. 
However, display values captured in the calibration process typically have 0-255 gray 
level values. 

A gamma value can be chosen as a standard target. For example, a gamma value 
of 2.2 is often assumed for uncorrected CRT displays and is reasonably linear with 
respect to L*, i.e., in the L*a*b* color space. As described in greater detail below, 
calibration module 24 may implement techniques that can match this gamma behavior 
approximately to within an expected error. 

The matrix of EQUATION 1 can be redefined using standard values of 
chromaticity for R, G and B and for Yxy (luminance and chromaticity) of the white point 
of CRT 12. The matrix M for converting RGB to XYZ can be defined according to the 
chromaticities of the RGB values as follows: 



EQUATION 3 
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In order to achieve a particular targeted luminance for white Y WP} in units of 
candelas/meter 2 , the matrix M can be multiplied by Y wp . Thus, EQUATION 1 can be 
rewritten as: 
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where the values of RGB chromaticities are assumed to be constants. 

Experimental data gathered on several dozen CRTs, ranging in age from zero to 
eight years old, has indicated that the RGB chromaticities tend to remain very constant 
over time and between brands of CRTs. Thus, the variables that require adjustment are 
primarily the Yxy white point and the RGB gamma curves of the CRT. However, more 
complex adjustments which include changes in RGB chromaticities can also be 
performed as desired. 

The measured behavior of an uncalibrated CRT can be defined as follows: 
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EQUATION 8C 



The vector XYZdcis the dark current offset introduced by the measurement device 
plus stray light. Typical values for the X-Rite DTP92 measurement device are on the 
5 order of 0.3 candelas/meter 2 , compared with the white point luminosity of 8 1 .5 
candelas/meter . 

The values of RGB02 represent the black onsets, which indicate the 
threshold values of RGB above which an increase in XYZ can be detected. For 
perfect CRT behavior (as represented in EQUATION 1) the values of RGBo2are zero. 
1 0 The vector XYZdc may by a constant or may vary as a function of (1 .0 - a) where: 

less 

£ EQUATION 9 

R a = (XYZ (RGB)- XYZ(RGB=0))/(XYZ(RGB=255) - XYZ(RGB=0)) 

m 

(ft 15 For example, for some measurement devices, the value of XYZ for RGB=0 can vary 

significantly, while the value of XYZ for RGB=255 is very stable and repeatable. Thus, 
Jf 5 for some measurement devices the value of XYZdc can be multiplied by (1 .0 - a). 
W FIG. 3 is a flow diagram illustrating a calibration technique according to an 

f=3 embodiment of the invention. As shown, calibration module 24 initializes LUT 29 within 
^ 20 video card 26 (31). In other words, LUT 29 is set to linear in order to ensure that the 
output of CRT 12 is in a known state. Exemplary pseudo code for the initialization is 
provided below. 



^define numEntries 1024 
25 Mefine maxl6bit Oxffff 

unsigned short videoLUT [numEntries] ; 
int i; 



(for i=0; i< numEntries ;i++) 
30 videoLUTfiJ = (unsigned short) (0.5 + maxl6bit * ((float)i/(float) (numEntries- 1)); 
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Although this linear video LUT will be altered during the calibration procedure, 
it is advantageous to ensure that videoLUT[0] is always kept at 0. This will ensure that 
the measured value of RGB=0 is a constant throughout the calibration process. 
5 After the LUT has been initialized (3 1) 5 the settings of CRT 12 can be adjusted 

(32), either automatically by calibration module 24 or manually by a user. Adjusting the 
analog settings of CRT 12 further ensures that the initial state of CRT 12 is known and 
defined. As mentioned above, the state of the RGB target space can be defined by white 
point and gamma values for RGB, assuming the RGB chromaticities are constant. The 
10 values of the analog settings are typically referred to by the following names. Contrast 
refers to the setting that adjusts primarily the luminance of the near white colors for the 
CRT (which affect RGB values together), and secondarily affect dark luminance values. 
Brightness refers to the setting that adjusts primarily the luminance of the dark luminance 
colors for the CRT (which affect RGB values together), and secondarily affects near 
15 white values. Gain refers to the setting that adjusts similar values as contrast, but does so 
separately for each R, G, and B channel. Bias refers to the setting that adjusts similar 
J3 values as brightness, but does so separately for each R, G, and B channel. 

The above analog settings can be set either manually, e.g., using input buttons, 
dials and the like presented by CRT 12, or automatically via direct communication 
M. 20 between computer 14 and CRT 12. An example of such direct communication is the 
Video Electronics Standards Association (VESA®) standard known as ddc/ci. In one 
specific implementation, the values of contrast, brightness, and the gains and biases of 
R,G, and B can be adjusted to achieve the following measured behavior: 

1) The value of Yxy for white (RGB=255) is set as close as possible to the 
25 target values of Yxy with the exception that the value of Y is set slightly (i.e. 1 %) 

higher than the value of Y for the white of the RGB target space. This allows 
room to perform remaining corrections via LUT 29 in video card 26. The target 
value of Yxy for white can be chosen by the user according to the preferences or 
needs of the user. 
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2) The value of x,y for a dark gray (e.g. RGB=50) is set as close as 
possible to the value of x 5 y for the target white. Again, how the target is defined 
may depend on the preferences of the user. 

3) The luminosity Y of the dark gray can be adjusted in one of two ways 
5 that are essentially equivalent. In one case, the dark gray luminosity is set such 

that the value of brightness is as high as possible while still measuring perfect 
black at RGB=0. "Perfect black" is defined as the value of Y for the CRT with 
RGB=0, and all analog adjustment values set to minimum. In the other case, the 
highest value of brightness is determined which results in a value of the black 
1 0 onsets of RGB02 = 0 or slightly above 0. 

For some CRT's, the gamma value is approximately 2.35. In that case, it has 
been found that the desired value of brightness occurs when L* = 13.0 - 14.0 at RGB=48. 
The relative values of RGB bias can be adjusted to measure a*=b*=0. In that case, all 
calculations for CIELAB assume that XnYnZn is determined by the value Yxy of the 
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% 1 5 white point of the RGB target for CRT 1 2. 

2 After adjusting the settings (32), output of CRT 12 is measured (33). For 

yi 

y3 example, calibration module 24 may direct the user to attach a calibration device to the 
i screen of the CRT to measure output. When the measurement device is attached, 

r — 

nJ calibration module 24 can perform a calibration routine that involves displaying and 

i 3 I 

y= 20 measuring color samples. For example, data capturing may involve displaying and 
F measuring RGB gray values ranging from 0 to 255. All of the gray levels can be 

measured, or alternatively, a subset of gray levels can be measured in order to reduce the 
time it takes to capture data. Measurements in 1 5 gray level increments can provide a 
suitable level of accuracy. Each gray level defines a unique "neutral color." Neutral 
25 colors refer to colors having substantially equivalent color values. For example, neutral 
colors for a CRT are colors having substantially equivalent values for the R, G, and B 
channels. 

The calibration routine may display a circle representing RGB neutral color 
(typically a circle with a ten centimeter diameter or larger) against a black (RGB=0) 
30 background. The black background can help improve the ability to measure the darker 
colors. A number of XYZ values can be captured for each RGB value that is displayed in 
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order to obtain an estimate of measurement variability, i.e., an error value (sigma) 
associated with the measurement device. The calibration routine can then proceed to 
display the next RGB value and repeat the process until all RGB values, or a subset 
thereof, have been captured. 

One concern in measurement of dark emissive colors is that some 
measurement devices clip the data. In other words, for some measurement devices, 
values of X, Y, or Z in the range of RGB = 0 to RGB=40 may be returned as 0 when a 
truly accurate value would be finite but small, e.g. 0.4 candelas/meters . To compensate 
for clipping in the measurement device, it may be advantageous to display a substantially 
white trace during the calibration routine to slightly bias the measurement device. The 
white trace may be a fixed white circle (or halo) having color values of RGB=255. The 
scattered light from this halo is a constant non-zero quantity and can offer the beneficial 
effect of adding a slight bias to the values of XYZ measured in darker colors. Hence, all 
values of RGB from 0 to 255 measured in this manner can be guaranteed to be within the 
dynamic range capabilities of the measurement device. 

FIG. 4 illustrates one example of an imaging station 10 displaying a white trace 
during calibration. Measurement device 42 can be affixed to the display screen of CRT 
12, and substantially white trace 44 can be displayed during the calibration procedure to 
bias the measurements in order to ensure that the values of RGB from 0 to 255 are within 
the dynamic range capabilities of measurement device 42. The color being measured 
may fill the interior of white trace 44, or alternatively both the color being measured and 
a black background may be included in the interior of white trace 44. In either case, 
white trace 44 can be used to provide a desired amount of biasing. For example, biasing 
the output measurements by approximately 0.5 to 1.0 percent is usually sufficient. 
However, depending on the measurement device, a larger amount of bias may be useful. 
White trace 44 is illustrated as having a halo shape, but could be implemented with any 
shape or design. Also, biasing with other colors may be useful for some measurement 
devices. 

After the output has been properly measured (33), the parameter values of the 
device model are chosen (34). The number of adjustable parameters can be defined in the 
device model to be less than a number of measurements used to define the measured 
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output of cathode ray tube 12. By defining the number of parameters to be less than, or 
substantially less then the number of measurements used to define the measured output, it 
can be ensured that the calibration process will not overcompensate for measured results. 
In other words, there is always some level of error that can be attributed to the 
measurement device or other external factors. The invention can avoid improperly 
compensating for external factors unrelated to the actual device output by limiting the 
number of adjustable parameters to be less than, or substantially less then the number of 
measurements used to define the measured output. Specifically, the invention can avoid 
adjusting each input according to the measured output for that given input. Rather, the 
invention may implement a device model based partially in theory and partially on 
chosen parameter values. In this manner, adjustments do not overcompensate for 
measured errors that should be expected. Values for the adjustable parameters can be 
chosen during the calibration process in a manner that minimizes the error between 
analytical expected outputs and measured outputs. 

Which parameters are defined as adjustable may depend on the implementation. 
Indeed, for different imaging devices, the adjustable parameters would typically be 
different. However, in the example of calibrating CRT 12, the adjustable parameters may 
be a gamma value and a black onset value. The gamma value refers to the value of a well 
know parameter in the art (y) that indicates the rate of change in light intensity with 
change in digital device value. The black onset value is another well know term in the 
art, which refers to the point at which there is a measurable increase in light. 

The parameter values of the device model can be chosen by an error minimization 
process in order to define a robust device model that is accurate to an extent approaching 
theoretical limits. In one example, the device model of a CRT is defined using the 
following process. First, the captured XYZ data is adjusted with respect to a dark current 
bias. In particular, the dark current bias (XYZ dc defined above) of the XYZ data can be 
subtracted from the raw measured XYZ data. If the value of XYZ dc is a constant 
for a particular measurement device or method of measurement, then the value of XYZdc 
can simply be subtracted from all values of XYZ. However, if the value of XYZdc varies 
according to the XYZ value, the adjustment to the XYZ data should be performed 
accordingly. The result of this adjustment to the XYZ data is that the value of XYZdc in 
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EQUATION 7 is zero. Thus, only the remaining parameters must be determined in order 
to accurately characterize the CRT behavior. 

Next, the white point values are set. For example, the values of (Yxy) wp 2 can be 
plugged into EQUATION 7 above. In that case, for (Yxy) wp 2, the values of XYZ wp 2 are 
5 given by the XYZ data at RGB=255, and the definitions for x and y are: 

x =X/(X+Y+Z) EQUATION 10A 

y = Y/(X+Y+Z) EQUA TION 10B 

10 Next, an error function can be defined that depends on the chosen parameter 

values, in this case the gamma and black onset parameter values. For example, an error 
function can be defined as follows: 



EQUATION 12 



Error (y r2 ,y g2 ,y h2 , RGB) = £ AE? (y r2 , y g2 , y b2 , RGB o2 ) 
J =° EQUATION 11 

5 15 

si 

yj where 
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h AE?(y r2 ,y g2 ,y b2 ,RGB o2 ) = (I * (RGB/ , y r2 , y g2 , y b2 , RGB o2 )-L*,) 2 + 

M 20 (a^RGB^y^y^y^RGB^-a*,) 2 + 

(b * (RGB, ,y r2 ,y s2 ,y h2 , RGB o2 )-b* i ) 2 



L*(), a*() and b*() are the analytical expected values of L*a*b*, whereas the values 
L*i , a*i , b*i are the measured values of L*a*b*. The analytical expected values and 
25 measured CIELAB values can be derived from analytical and measured XYZ. The 

analytical expression for X,Y, and Z, as a function of RGB, XYZ dc , y t 2, Yg2> Yb2, RGB 0 2, is 
defined in EQUATION 6 and EQUATION 7 above. The calculation for CIELAB 
values from XYZ is given by: 
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I* = 116/(r/yj-16 EQUATION 13A 

a* = 5Q0[AX/X n )-f(Y/Y„)] EQUATION 13B 

-6* = 200[/(7/yj-/(Z/ZJ] ^IMIWW 7JC 



where 



/(") = 



{w) m for w> 0.008856 ^ 1 
7.787(w) + 16/116/orw<0.008856 J EQUATION 14 



1 0 The values of RGBi and X YZi are the corresponding values of displayed RGB and 
measured XYZ for each input value of RGB ranging from 0 to 255. The summation 
performed above from i=0 to i=N-l assumes a list or array of (N) XYZ measured values 
and their corresponding N values of RGB used to generate the measured color. RGB=0 
may be excluded in this calculation due to the possibility of a singularity or non-linearity 
*0 15 of behavior in the region near or at RGB=0. In other words, not including RGB=0 in this 

3 

u calculation can ensure that the overall RGB to XYZ behavior of the CRT behavior is 

jjfi being characterized without being skewed by behavior at or near RGB=0. 

H 8 At this point, values for the parameters can be determined by error minimization. 

o 

For example, an error minimization can be performed on the adjustable parameter using 
20 any suitable error minimization method such as a chi-squared method or a least squares 
fit. Other error minimization techniques could also be used. In any case, nominal initial 
values for the chosen parameters can be used (in this case gamma = 2.2 and black onset = 
0.05). After the error minimization has been performed, calibration module 24 may 
check the average and maximum error between the measured data and the analytical 
25 expected behavior that has been calculated. 

Experiments have shown that typical values for average error range from 0.3 - 0.5 
delta e and that maximum errors range from 0.6 - 1 .0 delta e. These errors appear to be 
due to systematic quantization errors in the analog-to-digital (A/D) and digital-to-analog 
(D/A) circuitry of the video card and the measurement device causing +/- shifts in XYZ 
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as a smooth function of RGB. The sigma of noise of measurement for a particular RGB 
color is typically quite low, e.g., 0.1 delta e. Thus, it should be expected that this 
quantization error should be evenly distributed about the perfect gamma curve behavior. 
Hence, although the average and maximum errors from analytical expected behavior may 
be slightly larger than variability due to random measurement noise, one should expect 
that the sum of the residuals should be close to zero. 

To even further improve the error minimization, additional techniques can be 
employed. For example, error minimization routines sometimes do not find the best 
minimum for the error function. In particular, a fatal error can occur in the minimization 
routine, such as divergence rather than convergence. To avoid this, it may be 
advantageous to find the best estimates possible for the initialization of the parameter 
values in order to ensure that the error function is as near minimum as possible to its 
ultimate value. 

One way to pre-optimize the parameter values discussed above is to measure a 
smaller set of data for each channel separately and to calculate and minimize the error of 
the separate channels first. In this approach, a smaller set of single channel data is 
obtained, i.e., the other two channels are set to 0 during data capture. Additionally, the 
predictions for L*a*b* can be performed with two of the three RGB values set to 0. 
Finally, the error minimization can be performed only on the parameters for that one 
channel. These additional techniques can better insure that the error minimization 
technique will indeed identify parameter values that minimize the error between 
analytical calculations and measured behavior. 

Once the parameter values of the device model are chosen (34), calibration 
module 24 can create and load entries into LUT 29 (35). In this manner, LUT 29 can 
adjust the image data to account for drift in CRT 12. LUT 29 may be in the form: 



R} =fr ( R ) = R <>2 + (^max2 ~ R ol) R 



Yr\lYrl 



EQUATION ISA 



G' =f g -\G) = G o2 +(G max2 -G o2 )G 



EQUATION 15B 



B'=f-\B) = B o2+ (B max2 -B o2 )B 



EQUATION 15C 
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At this point all the values of EQUATIONS 1 5 A- 1 5C are known with the 
exception of RGB max 2. To obtain the values for RGB max2 , the expressions for the 
uncalibrated CRT indicated by EQUATIONS 7 and 8A-8C can be equated with the 
expression for the desired CRT behavior indicated by EQUATION 6 after substituting 
R'G'B' for RGB in the expressions for the uncalibrated CRT. In other words, one can 
attempt to predict how close to the desired CRT behavior the uncalibrated CRT will be 
after the adjustments stored in video LUT 29 are applied. Substituting these expressions 
for R'G'B' into EQUATIONS 7 and 8A-8C yields the following: 







fx \ 


Y 










K^dc ) 



+ Y wpl M(x wpl ,Y wp \r sl (G) 

f>2 (*) J 



EQUATION 16 



where 



fr2 (*) = 



f'„2 (B) = 





-R02) 


Yrl 


_ (1.0- 


-K2) . 




( G nmx 


-gJ 




_ (1.0 
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-b o2 ) 


Ybi 


(1.0- 


-B o2 ) . 
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B 



EQUATION 17 A 



EQUATION 17B 



EQUATION 17C 



Equating EQUATION 1 6 with EQUATION 6 for the ideal CRT behavior with a 
gamma of 2.2, with additional measured dark current bias added yields the following: 

EQUATION 18 



A dc 




7hW 




(y \ 


Y dc 


+ Y wpX M{x wpX ,Y wpU 


/..«*> 






\^dc ) 




Ui(*)J 




K Z dc) 



+ Y wp2 M(x wp2 ,Y wp2/ 



f' S 2(G) 



f b2 (B) 



which can be reduced to: 
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EQUATION 19 



/'«(*). 



= -^L M (x wp2 , y wp2 )M (x wpl , y wp[ ) 



1 wp2 



The product of the matrices is diagonal, and results in the following three 
equations: 
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Q 

W 

01 

cn 

M- 10 



7 7 
Y Y 

1 wp2* \vpR2 



/,.(*) 



f' S 2(G) = 



F 7 

wp 1 wpG 1 

wp 2 wpG 2 



Y Y 

wp I ivpfl 1 r ( n\ 

y y Jb\\ D ) 

1 wp 2 1 wpB 2 



EQUATION 20A 



EQUATION 20B 



EQUATION 20C 



RGB maX 2 can now be obtained by substituting EQUATIONS 2A-2C and 
EQUATIONS 17A-17C into EQUATIONS 20A-20C as follows: 



_ (1.0- 
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JV2 


Y Y 

wp\ wpR\ g /r] 

Y Y 


EQUATION 21 A 


( G max 

(1.0- 


-gJ 

~G o2 ) . 


G Xg] = 


Y Y 

x wp\*wpG\ Q/ g i 

Y Y 

1 wp2 l wpG2 


EQUATION 21 B 


( 5 max 

_ (1.0- 


-5 o2 ) J 


B rbl = 


Y Y 

wp\ wpB\ gy b{ 

Y Y 


EQUATION 21C 



1 5 and solving for R max , G max and B max as follows: 



^nax=^2+(l-0-^ 2 ) 



Y Y 

Y Y 

1 wp2 1 wpR2 



EQUATION 22 A 
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EQUATION 22B 



EQUATION 22C 

The values of EQUATIONS 22A-22C can then be substituted into EQUATIONS 15A- 
1 5C to define entries of LUT 29. The procedure can be repeated as desired, each time 
using a correction function calculated from EQUATIONS 15A-15C from a previous 
iteration. In other words, the adjustments may be implemented as a closed loop such that 
the adjustment values continue to converge with each iteration of the procedure. The 
procedure may be repeated each day, or may be repeated as desired. For example, it may 
be advantageous to repeat the procedure before a contract proof image is displayed. 
Experiments have shown that typically CRTs drift by only 1 .0 to 3.0 delta e over an 
approximately 24 hour period. Thus, for many applications one calibration per day is 
adequate. 

Additionally, intelligent process control may be implemented such as by 
calculating the adjustment to the most recent correction function, and performing a partial 
correction rather than a full correction. The degree of correction can be adjusted 
based on the noise of measurement versus the magnitude of error. Additionally, 
intelligent process control may be implemented by correcting the most recent correction 
function only when the average and/or max errors exceed a predetermined value. 

The target Yxy value for white can often be achieved to an acceptable level by 
automatic adjustment of the CRT analog parameters as described above. Optionally, 
calibration module 24 can adjust the gammas and onsets without modifying the 
maximum correction values. Likewise, the onsets may be minimized by optimization of 
the CRT analog parameters. In that case, the calibration method described above can be 
used to adjust only the gamma behavior of the RGB channels in order to ensure 
accurate gray balance. Even if adjustments are performed on each RGB channel 
for the black onset, gamma, and max value, minimizing the magnitude 



G max =G o2+ (\.0-G o2 ) 



Y Y 

wp\ wpG\ 
^wp2^wpG2 



gi 



B max =B o2+ Q-0-B o2 ) 



Y Y 

Y Y 

1 wp2 1 wpB2 
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of adjustment for black onset and max value is desirable when limited by an 8 bit video 
LUT in order to minimize the effects of quantization. 

If the choice of the analog settings is used to calibrate the system white point to a 
value close to the target white value (for example, no more than +/- 0.5 delta e error in 
5 L*,a*, or b* using the target white value for X n Y n Zn) there are at least two ways to 

characterize the gamma curve properties of the system. One way is to set the white point 
of the analytical model to the value of the measured white point of the CRT system after 
the analog selections are performed. This should result in a gamma curve for R, G, and B 
that is minimized in error with respect to a "relative white" reference target. In other 
10 words a target that is exactly 2.2 gamma for R,G, and B, and a white point equal to the 
current white of the CRT, i.e. close to but not exactly the target value of Yxy. 

Another way to characterize the gamma curve properties of the system is to set 
jz the white point of the analytical model to the desired target value for the white point of 
Q the CRT system. When this approach is used, the calculated model will always have a 

1 . 3 

2 15 small error relative to the measured data at RGB = white. However, the overall model 

may obtain a better fit to the measured data. This means that small errors in the white 

m 

y3 point will to some extent be accounted for in the overall analytical model for the CRT. 
L Consequently, the final gamma curve of the system should have a minimized error with 

P respect to the "absolute white" reference target, i.e. the reference target that has 2.2 

UJ 

M 20 gamma for R, G, and B as well as a white point that is that desired target white point of 
rf the system. 

Additional setting conditions for the CRT that can be checked before performing 
the above procedure may include the following. First, all values of the RGB reference 
space in the equations above should be less than or equal to the values of Y for the 

25 physical RGB device being calibrated. Second, all values of the RGB onsets for the 

physical RGB device being calibrated should be in a "correctable" range. For example, 
with current CRT technology, an estimate for what constitutes "correctable" is 5 < RGB02 
< 40 in units of gray levels. Third, all values of XYZdc should be positive to avoid 
clipping, and should be below a "reasonable" level. In that case, with current CRT 

30 technology, an estimate of what constitutes "reasonable" is 0.5 candellas/meter 2 
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The first setting condition simply means that the maximum values of RGB can 
always be reduced in order to hit the target, but cannot exceed 100% if the intensity of a 
channel is too low. The second setting condition can ensure that a good black (RGB=0) 
can always be achieved. The third setting condition can ensure that the measurement 
device, plus stray light, plus CRT minimum luminance values are reasonably low. The 
CRT white point should be set close to the vicinity of the targeted white point such that 
the maximum adjustment values of RGB for the RGB device being calibrated are in the 
range of 0.95 < RGBmax < 1.0. 

FIG. 5 is another flow diagram illustrating a more general overview of a 
calibration technique for a CRT according to the invention. As shown, the output of a 
CRT is measured (51). Then, parameter values of a device model of the CRT are chosen 
(52). Finally, image data is adjusted so that output of the CRT is more accurate (53). 

In particular, when the number of adjustable parameters of the CRT is established 
to be less than a number of measurements of the CRT, calibration can be improved. For 
example, rather than adjusting each input of the CRT according to a measured output, 
i.e., a one-to-one mapping of input and measurements, the invention establishes a device 
model that has fewer adjustable parameters than the number of measured outputs. In this 
manner, the technique does not overcompensate for measured results that have errors 
which would be expected. In other words, even after applying an adjustment to input at 
the video card, one should expect small errors due to factors unrelated to drift in the CRT 
such as, inaccuracies in the measuring device. The invention can achieve a balance 
between analytical modeling and empirical measurement so as not to overcompensate for 
measured errors that should be expected. 

By using the techniques described herein, improved calibration of CRTs can be 
achieved. Specifically, the techniques described herein can be used to calibrate a display 
such that an average color error is approximately less than (0.75 delta e) from an 
analytical expected color output and a maximum color error is approximately less than 
(1.5 delta e) from the analytical expected color output. More specifically, average color 
errors approximately between (0.3 delta e) and (0.75 delta e) from the analytical expected 
color output and maximum color errors approximately between (0.6 delta e) and (1.1 
delta e) from the analytical expected color output can be achieved. Even more 
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specifically, average color errors approximately between (0.3 delta e) and (0.4 delta e) 
from the analytical expected color output and the maximum color errors approximately 
between (0.6 delta e) and (0.8 delta e) from the analytical expected color output have 
been achieved. Average colors errors approaching (0.3 delta e) approach the theoretical 
5 limits that can be achieved with the current state of CRTs and measurement devices. As 
measurement devices and CRTs improve, however, the theoretical limits of calibration 
may also improve. In that case, the invention may achieve even better average errors and 
maximum color errors. 

FIG. 6 is another flow diagram generally illustrating a calibration technique that 
10 can be applied to any imaging device in accordance with the invention. As shown, the 
imaging device is characterized with a device model (61). In that case, an average error 
between an analytical expected value of the device model and a measured output of a 
subset of device values of the imaging device can be made on the order of an expected 
error. Furthermore, in some cases, a subset of device values substantially corresponding 
^15 to neutral colors may comprise the output. After characterizing the imaging device (61), 



01 the image rendering by the imaging device can be adjusted to achieve a target behavior 

(62). For example, adjusting the imaging rendering comprises adjusting image data, such 



as by using a LUT or by creating dynamic color profiles as described below with 



5 

s . 

fU reference to FIG. 7. 

y 

1^20 The device model used to characterize the imaging device may include one or 

more adjustable parameters. In that case, the technique may further include choosing 
adjustable parameters of the device model. Also, the number of chosen parameters may 
be less than a number of measurements used to define the measured output of the 
imaging device. For example, in the case where the imaging device is a CRT, the 

25 adjustable parameters may comprise a gamma value and a black onset value. However, 
in other examples, the parameters can be specified according to the imaging device. 

FIG. 7 illustrates an alternative embodiment of imaging station 10. In this case, 
the same calibration techniques can be used to calibrate CRT 12. However, rather than 
apply the adjustments using a LUT, this embodiment creates dynamic color profiles in 

30 color matching module (CMM) 72. In other words, after calibration module 74 has 

performed the calibration procedure, the adjustment to account for drift in CRT 12 can be 
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applied by generating dynamic color profiles in CMM 72. The color profiles are 
dynamic in the sense that they change with calibration measurements. Thus, as CRT 12 
drifts, the dynamic color profile may change. In other words, CMM 72 may receive the 
- calibration data from calibration module 74 and can incorporate the calibration 
information into a device profile specifically for CRT 12. Adjusted image data can then 
be fed through display driver 75 and video card 76, possibly without any additional 
adjustments. In this manner, the need for a LUT within video card 76 can be avoided. 

FIG. 8 illustrates a soft proofing system 80. Soft proofing system 80 may 
implement one or more aspects of the invention to realize accurate color rendering and 
color matching in a soft proofing process. Soft proofing system 80 includes an 
administrative computer 82. Administrative computer 82 can be thought of as a server 
computer for soft proofing system 82. Administrative computer 82 may serve up images 
to imaging stations 10A-10D (hereafter imaging stations 10). Color specialists at 
imaging stations 10 can inspect the images, and possibly provide feedback by marking or 
highlighting the images and retuning marked-up copies to administrative computer 82. 
Upon receiving feedback, an administrator may implement changes to the image using 
administrative computer 82. Once the administrator and the reviewers associated with 
imaging stations 10 reach agreement on the appearance of the color image, the image can 
be printed via a printing press or another high quality printing device. 

Administrative computer 82 may be directly coupled to the various imaging 
stations 10, possibly forming a local area network (LAN), or alternatively, the 
administrative computer 82 may be coupled to imaging stations 10 via a wide area 
network or a global network 84 such as the Internet. The various imaging stations 10 
may implement one or more of the calibration techniques described above in order to 
improve calibration and thus improve color accuracy. Indeed, improved calibration can 
directly impact the ability to achieve an effective soft proofing system 80 or another color 
intensive application. 

In one embodiment, imaging stations 10 collectively define a set of CRTs, with 
each CRT in the set being calibrated to within the errors outlined above. In that case, it 
can be ensured that drift in the CRT of any imaging station 10 is adequately accounted 
for, such that contract proof color quality can be rendered at each imaging station 10. 
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Many aspects of the invention have been described as being at least partially 
implemented in software. Alternatively, exemplary hardware implementations may 
include implementations within a DSP, an application specific integrated circuit (ASIC), 
a field programmable gate array (FPGA), a programmable logic device, specifically 
designed hardware components, or any combination thereof. 

Although many aspects of the invention have been described in the context of a 
calibration technique for calibrating a display device in the form of a cathode ray tube, 
aspects of the invention may be readily applicable to calibration of other imaging devices, 
including liquid crystal displays, plasma displays, various printing devices, or other 
imaging devices such as scanner or the like. Accordingly, other implementations and 
embodiments are within the scope of the following claims. 
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